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Abstract 

We calculate the leading logarithmic QCD corrections to the decay Bg — > 77 in the 
two Higgs doublet model (2HDM) including O7 type long distance effects and estimate 
the restrictions of the 2HDM parameters, tanf3 and mn, using the experimental data of 
B — > Xg^ decay provided by the CLEO Collaboration. A lower bound for the charged 
Higgs mass rriH as a function of the renormalization scale fi is given for 2HDM model 
H. We further present the dependencies of the branching ratio Br{Bs — > 77) and the 
ratio l^+p/l^-p on niH and tan/3 including leading logarithmic QCD corrections. The 
dependence on the renormalization scale is found to be strong for both ratios. An addi- 
tional uncertainty arises from the variation of the parameters of the bound state model, 
{nib, ^s)- We see, that to look for charged Higgs effects the measurement of the branching 
ratio Br{Bs 77) is promising. 
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1 Introduction 



The experimental discovery of the inclusive and exclusive B — > Xg'y ^ and B —>■ K*'~f [0 
decays stimulated the study of rare B meson decays as a new force. These decays take place via 
flavor- changing neutral current (FCNC) b ^ s transitions, which are absent in the Standard 
Model (SM) at tree level and appear only through loops. Therefore, the study of rare B-decays 
can provide a sensitive test of the structure of the SM at loop level and may shed light on 
the Kobayashi-Cabibbo-Maskawa (CKM) matrix elements and the leptonic decay constants of 
the B-mesons. At the same time, such decays are in a very promising class to search for new 
physics beyond the SM, like two Higgs doublet model (2HDM), minimal supersymmetric model 
(MSSM), etc. 1^. Currently, the main interest is focused on these decays for which the SM 
predicts large branching ratio and which can be measured in near future in the constructed 
B-factories. The Bs '~f'~f decay belongs to this category. In the SM, the branching ratio of 
Bs 77 decay is of order 10^^ without QCD corrections. Including leading log (LLog) QCD 
corrections, the branching ratio {Br) b —>■ S77 is of the same order of mag nitude ~ lO'^ like 
b^ sl+l~ [|. 

The investigation of B^ 77 decay is interesting for the following reasons: 

• It is well-known, that the QCD corrections to the b ^ s'~f decay are considerably large (see 
[^] - 1^ and references therein). Therefore, one can naturally expect that the situation is 
the same for the b 577 decay. Recently the QCD corrections in the LLog approximation 



to this decay have been calculated and found to be large as expected - ||rT| . Note, that 



in the literature this decay without QCD corrections was analysed in the SM |T2[ - ||13 



and in the 2HDM P 



In Bg 77 decay, the final photons can be in a CP-odd or a CP-even state. Therefore 
this decay allows us to study CP violating effects. 

From the experimental point of view, Bg 77 decay can be easily identified by putting 
a cut for the energy of the final photons, e.g., the energy of each photon is larger than 



100 MeV . In this case, two hard photons will be easily detected in the experiments [|T^ . 
Finally, this decay is also sensitive to the physics beyond the SM. 



In an earlier analysis [|I^, the Br{Bs 77) in the 2HDM without QCD corrections was found 
to be enhanced with respect to the SM one for some values of the parameter space. In the 
present work, we study Bg 77 decay in the 2HDM with perturbative QCD corrections in 
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LLog approximation. In contrast to [0, who used the constituent quark model, we impose a 
model based on heavy quark effective theory for the bound state of the meson Bs- Further we 
perform an additional analysis with the inclusion of long-distance effects through the transition 
Bs ^ (pj ^ 77, which we call Oj-tjpe throughout this paper, see for details. We find, 
that the theoretical analysis is shadowed by large uncertainties due to the renormalization 
scale fi and the parameters of the bound state. The decay Bs 77 is dominated by the 
Wilson coefficient Cj-^'^ (see section 2), which is restricted in our analysis by the B Xs'j 
branching ratio provided by CLEO data Br{B Xs'~)) oc \C^^^\^^ see section 3. Without 
any improvement from the theoretical side, we see that the only chance to detect a deviation 
from the SM in Bs — * 77 decay lies in a possible enhanced branching ratio, which can be at 
most 1.4 ■ 10-6 in the SM and 2.1 ■ 10"^ in model II of the 2HDM (for mn = 480 GeV and 
large tanf3) resulting from our analysis, at ^ = 2.5 GeV including the Oy-type long distance 
effects. 

The paper is organized as follows: In Section 2, we give the LLog QCD corrected Hamil- 
tonian responsible for the b S77 decay. We further calculate the CP-odd A~ and CP-even 
amplitudes in an approach based on heavy quark effective theory, taking the LLog QCD 
corrections into account. In Section 3, we study the constraint analysis for the 2IIDM parame- 
ters rriH and tan(3, using the measured data on the branching ratio of the B Xg'y decay [Q]. 
Section 4 is devoted to an analysis of the dependence of the ratio [A'^'p/IA-p and the Br on 
the parameters fi (scale parameter), tan(3 and rriH and our conclusions. 

2 Leading logarithmic improved short-distance contri- 
butions in the 2HDM for the decay Bg 77 

Before discussing the LLog QCD corrections to the Bg 77 decay, we would like to remind 
the main features of the models which we use in further discussions. In the current literature, 
mainly two types of 2HDM are discussed. In the so-called model I, the up and down quarks 
get a mass via the vacuum expectation value (v.e.v) of only one Higgs field. In model II, the 
up and down quarks get mass via v.e.v of the Higgs fields Hi and H2, respectively, where Hi 
{H2) corresponds to first (second) Higgs doublet of the 2HDM. Note, that in this sense the 
Higgs sector of model II coincides with the MSSM extension of the SM. In the 2HDM, there 
exist five physical Higgs fields, namely, two charged H^ and three neutral Higgs bosons. The 
interaction Lagrangian of the quarks with the charged fields, which we need for the calculation 
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of the b S77 decay amplitude, is HT? 



^ ^ y ~ rndX'uiLdj\VijH^ + h.c. , (1) 

where L and i? denote chiral projections L{R) = 1/2(1 =1= 75) and ^ and ^' are the ratios of the 
two vacuum expectation values, vi and V2 of the Higgs fields Hi and H2, respectively. Vij are 
the elements of the CKM matrix. In model II, 

C = -1/^ = -tanP = -V1/V2 , (2) 

and in model I 

e = ^ = cotp = V2/vi. (3) 

After this preliminary remark, let us discuss the LLog QCD corrections to the b sj'-f 
decay amplitude in the 2HDM. The framework to incorporate short-distance QCD corrections 
in a systematic way is that of an effective low energy theory with five quarks, namely u,d,s,c,b 
quarks. The effective Hamiltonian is obtained by integrating out the heavier degrees of freedom, 
i.e. the top quark, and bosons. In the effective theory, only the lowest (mass) dimension 
operators, which are constructed by quark and gauge fields, are taken into account, since higher 
dimensional operators are suppressed by factors 0{ml/m^) and 0(m^/m^). 

The LLog QCD corrections are done through matching the full theory with the effective 
theory at the high scale /z = m\y and then evaluating the Wilson coefficients from down 
to the lower scale fi ~ 0{mh). In this way the LLog QCD corrections for the b —>■ s'yj decay in 
the SM are calculated in ^ - ||TT |. 



The effective Hamiltonian relevant for our process is 

-Heff = -^%v,,v,: jz c,{^m^) , (4) 

V2 i=i 

where the Oi are operators given in eq. (^ and the Ci are Wilson coefficients renormalized at 
the scale /i. The coefficients can be calculated perturbatively and the hadronic matrix elements 
< \^|Oi|i? > can be calculated using some non-perturbative methods. 
The operator basis of Ti-eff is given as 

01 = {sLal,,bLa){cLpl^CLp), 

02 = {sLalfibLp){cLpl''CL 
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O3 = (sLal^bLa) J2 (QL/bY qLls) 
q=u,d,s,c,b 



q=u,d,s,c,b 

O5 = (sLalfibla) {QRPI^QRP 
q=u,d,s,c,b 

q=u,d,s,c,b 

e 



^7 = T^^a(^^^'^(^bR + msL)baJ^^'', 



Os = -^s^npa,,im,R + msL)bpg'''''', (5) 

where a and (3 are SU{3) colour indices and JF^*^ and Q'^'^ are the field strength tensors of the 
electromagnetic and strong interactions, respectively. 

For the reason given below, the LLog QCD corrections can be calculated in analog to the 
SM. In the 2HDM, the charged Higgs fields are present and give new contributions due to the 
their exchange diagrams. Since the interaction vertices of the charged Higgs bosons and quarks 
are proportional to the ratio mq/mw, where rriq is the mass of the quark and mw is the mass of 
the W boson, the main contribution comes from the interaction with the t-quark. We neglect 
the contributions coming from u and c quarks, since their masses are negligibly small compared 
to mw- In this case the calculations show, that the new contributions modify only the Wilson 
coefficients C7 and Cg of the operators O7 and Og at mw scale and do not bring any new 



operators |Tj]. Therefore the operator basis used in the 2HDM is the same as the basis used 
in the SM for the b — > 577 decay. 

Denoting the Wilson coefficients for the SM with Cf^^(mpi/) and the additional charged 
Higgs contribution with C^ijny/), we have the initial values p!8| : 

CflAmw) = 0, 
C^'\mw) = 1 , 

^o.r, , 3x^-2x2 _8a;3-5x2 + 7x 



4(x - 1)4 24(x - 1; 
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/■ \ 3x X -\- 5x -|- 2x , , 



and ig, 1T§ 



Ctsimw) = 0, 



6(1/ -1)3 - 12(y-l)2^ ' 
+ ^tl2(^^^^+ 72(,-l)3 ] 
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2(^-1)3 - Aiy-ir 

A{y-iY ""y^ 2% -1)3 



where x = m'^/m'^ and ?/ = m^/m\j. Here the parameters ^ and ^' are given in eqs. (0) and 
(§). From eqs. and (0) the initial values of the coefficients for the 2HDM are defined as: 

Cr^'^'imw) = 1, 

C^^'^'^mw) = Cf^(mvK) + Cf(m^) , 

t^8'™'(^H^) = Ci^'{mw) + C^imw). (8) 

Using the initial values of the Wilson coefficients Cf^^^ , we can calculate their contribu- 
tions at any lower scale as in the SM case. Here we would like to make the following remark: 
Since in our case there exists a charged Higgs boson with a mass larger than mw, the correct 
procedure to calculate the Wilson coefficients at a lower scale /i has two stages: First, we cal- 
culate the value at mw starting from mn and second, we evaluate the result from mw down to 
a lower scale fi. We assume that the evaluation from m^ to mw gives a negligible contribution 
to the Wilson coefficients and therefore we consider only their evaluation between mw and a 
lower scale fi. 

Using the effective Hamiltonian in eq. M), the amplitude for the decay Bs 77 can be 



written as [|] - 

A{Bs ^ 77) = A+T^^T^" + lA-T^.f^'' , (9) 

where T^.^ = |e^^„/3^"^. Here {A') is the CP -even (CP-odd) part in a HQET inspired 
approach : 

^+ ^ aernGp ^ A m|;^ (mf ^ - ^) ^^^^ 
V27C m\ *1^3A,(mB. -A.)K-^^ + m^-^0 ' 
■{-mbJ{mb) +msJ{ms))D{n) , 



9 ml^^ + 



V27r \3mBAs{mB, - As) 

- Y.Qpimg)Cgifi)+ KA(mb)+m,A(m,))D(/i)) , (10) 

where Qq = | for g = m, c and Qg = — | for q = d,s,b. In the calculations, we have used 
the unitarity of the CKM-matrix J2i=uct ^^^6 = and have neglected the contribution due to 



KTs^wb ^ ^t*s^tb = At. The function g_ is defined as: 

The parameter enters eq. ( p!OD through the bound state kinematics p. ml^^ and m^-^-^ are 
the effective masses of the quarks in the Bs meson bound state , 

{nfj^f = p'^ = ml- 3A2 , 

(mtfff = p'^ = {mlfff-ml+2msKs, (12) 



where A2 comes from the matrix element of the heavy quark expansion |19|. The LLog QCD- 



corrected Wilson coefficients Ci.,,6(At) - fTll enter the amplitudes in the combinations 



Cu{^^) = Cdifi) = {CM - CM)Nc + C4(/i) - C6(/i) , 

c,(^) = (Ci(^) + C3(/i) - c,{ii))K + CM + CM) - CM , 
CM = CM = {CM + cMKNc + 1) - n,cm) - cm , 

= CM + CMNc, (13) 

where Nc is the number of colours {Nc = 3 for QCD). While Ci...6(/u) are the coefficients of the 
operators Oi,,s, C^^^ {fi) is the "effective" coefficient of O7 and contains renormalization scheme 
dependent contributions from the four-quark operators Oi,,s in ^e// to the effective vertex in 
b 57. In the NDR scheme, which we use here, Cj^^ {^) = Ct{^) — ^C^{^) — Cq{^)^ see [|18| 
for details. The functions I{mq), J{mg) and A{mg) come from the irreducible diagrams with 
an internal q type quark propagating and are defined as 

2 

771 

I{77ig) = 1 + ^A(m,) , 
m 



Bs 

771%, — 4m^ 
J(m,) = l-^-^A(m,), 





Am'i 

) 






J4:77ll- 
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A{7r^g) = |ln( ^ ^ ) - ^vr for ^ > 1, 



A(m„) = - 2arctan(-!^ -) - tt for — ^ < 1. (14) 

77iBs I 4m^ 

In our numerical analysis we used the input values given in Table (|l|). 
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-1 
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At 


U.U4 


F (' R ^ 


4.uy ■ iU (^LreV ) 


f 

jBs 


u.z [Kje V j 


77?, D 

'"'tis 


5 369 fGeV) 


TD ± 
II it 


17^ mpV"! 

X 1 1 V.J V V ( 




80.26 (GeV) 


mz 


91.19 (GeV) 


a(5) 


0.214 (GeV) 


as{mz) 


0.117 


A2 


0.12 (GeV^) 



Table 1: Values of the input parameters used in the numerical calculations unless otherwise 
specified. 

3 Constraint analysis 



There is a considerable interest in the constraints of the parameter space of the 2HDM, espe- 
cially in model II, since its Higgs sector coincides with the minimal supersymmetric extension 
of the SM one. The free parameters of the 2HDM are the masses of the charged and neutral 
Higgs bosons and the ratio of the v.e.v. of the two Higgs fields, denoted by tan j3. In our 
analysis the neutral Higgs bosons are irrelevant, since they do not give any contribution to the 
b — s> S77 process. Therefore we consider as free parameters the mass mn of the charged Higgs 
boson and tan(3. By using existing experimental data, it is possible to find restrictions on the 
parameters mn and tanjS. 

The model independent lower bound of the mass of the charged Higgs mn > 44 GeV comes 
from the non-observation of charged H pairs in Z decays ||2^ . There are no experimental upper 
bounds for tuh except mn < 1 TeV to satisfy the unitarity condition ||21|. For model II, the 
constraints have already been studied. Top decays give mu > 147 GeV for large tan/5 



The lower bound of tan(3 is found to be 0.7 due to the decay Z ^ bb ||23| and in addition for 
large tan(3 the ratio tan(3/mH is restricted. The current limits are tanf3/mH < 0.38 GeV^^ 
1 2^ and tanf3/mH < 0.46 GeV^^ |^^, which come from the experimental results of branching 
ratios of the decays B —>■ ru and B Xtu. Recently, the exclusive decay mode B Dtv 
has been studied for model II. It was shown that this decay could be used to put an upper 
bound on tan(3/mH, with the sufficient data and the reduction of theoretical uncertainities. 
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Under these conditions, the upper bound was estimated as tanP/niH = 0.06 GeV~^. 

In the present work, we estimate the constraints for the 2HDM parameters using the result 
coming from the measurement of the decay B Xg'-f by the CLEO collaboration 

Br{B X,7) = (2.32 ± 0.57 ± 0.35) ■ 10"^ (15) 

To reduce the b-quark mass dependence let us consider the ratio 

BrjB ^ X,7) 
Br{B ^ X,eu,) 



R 



|Vcb|2 TTg{z) 

where g{z) is the phase space factor in semileptonic b-decay, 

g{z) = 1 - + 8^^ + / - 24:zHn z (17) 

and z = mc/rrih- 

Now we want to discuss the theoretical uncertainties present in the prediction of R. 

\v*v P 

• The ratio of the CKM matrix elements I'^^p = 0.95 ± 0.04 has an uncertainty which 
comes from the CP violating parameter [^]. 

• The function g{z) has an uncertainty coming from the masses mi, and mc via the ratio 
z = rric/mb. HQET provides a mass relation 



mb-mc = {niB - m^) [1 - + 0(^)] , 

2mBiTiD rn'l 



where Q = b and c , and rrif) are spin averaged meson masses, = 5.31 GeV and 
mfj = 1.97 GeV [^. Here Ai is the non-perturbative parameter, which characterizes the 
average kinetic energy of the b-quark in B meson and its value is obtained by QCD sum 
rules. Using the theoretical estimate for Ai = —(0.4 ±0.2) GeV^ the mass difference 



and the error quoted are given as 

- m, = (3.40 ± 0.03 ± 0.03) GeV. (19) 

Here the first error is due to the uncertainty in Ai and the second one is from higher order 
corrections. We take the central value of b-quark mass as = 4.8 GeV. The uncertainty 
in rrib is Amf, = ±0.1 GeV. Using the HQET result, we estimate the uncertainty in rric as 
Arric = ±0.16 and we get the error in g{z) as Ag{z) = 0.096 and Ag{z)/ g{z) = 17.8 %, 
where z = mc/mb = 0.29 is the central value. 
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The Br for the semileptonic B X^eVf. is 



Br{B ^ X^ev^) =Q.lQ2,±Qm. (20) 

Both the theoretical uncertainties and the experimental errors, as given in eq. (|T^) and 
result in an uncertainty in Cj^^ . Using 

Sw,(S^X,7) = 3.24-10-^ , 
Brmin{B ^ Xsi) = 1.40 ■ 10-^ , 

we get a possible range for \Cy^^\ as 

0.1930 < \Cj^^\ < 0.4049 . 

In the SM and 2HDM model II is C^^^ < 0, but in possible extensions it can be positive. 

Now we present the lower bounds of rriH for different values of the scale fi in Table (|^) for 
model II. We restricted \Cy^^''^^^^\^)\ by using the limits given in eq. (|2^). For model I, a 



(21) 



(22) 



mHmin [GeV] 


/i [GeV] 


480 


2.5 


302 


5 


235 


10 


158 





Table 2: The lower bounds of the Higgs mass ttih for different scales /i in model II. 

lower bound for the Higgs boson mass is absent. 

The parameter tanP has bounds strongly depending on the scale /i and rriH- In fig. (|l|), 
we plot the parameter tanjS with respect to rriH for 3 different fi scales (2.5, 5, 10) GeV in 
model II, by fixing \Cj-^^''^^^^^\ = 0.4049. We see, that the dependence of tanP {ran) on ran 
{tan(3) becomes weak for large values of itih {tan(3) and that a decreasing /i scale causes the 
allowed region in the tan(3 - mn plane is to be small. It is interesting that at yU = 2.5 GeV, the 
solution for tan/S - mu exists only in the region 0.4047 < \Cj^f''^"^''"^\ < 0.4049. Therefore the 
solid curve in Fig (|I]) is almost the allowed region for {tan(3,mH). For /i = 2.5 GeV, we get an 
empirical expression for the restricted region of the parameter set {tanjS, mn), 

C2 



with Ci = 
(e ^ 1 ) and mn 



tanP = Ci + 
0.067, C2 = 6.9 GeV^/"^ and m„ 



mm 



rrir 



(23) 



e. Here e is a positive small number 



480 GeV. 
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In the following analysis we restrict the coefficient \Cj^'^''^^^^''\ in the given region and study 
the resulting mn, tanf3 and scale /x dependencies of the ratio and the Br for the 

decay Bg — 77. 

4 Discussion 



In the rest frame of the Bg meson, the CP = —1 amplitude A~ is proportional to the perpen- 
dicular spin polarization el x €2, and the CP = 1 amplitude A'^ is proportional to the parallel 
spin polarization el. €2- The ratio is informative to search for CP violating effects 

in Bs —>■ 77 decays and it has been studied before in the literature in the framework of the 



2HDM without QCD corrections |T^. In our analysis we use three sets of parameters (m;,, Ag) 
given in (Table (^), which model the bound state 0. However, we do not present the figures 
for the ffist two. Here we analyze the LLog /i, and 2HDM parameters {rriH, tanjS) dependence 
of the ratio |y4"''p/|y4~p and present the results in a series of graphs (fig. 2-7). 

In fig. d^) and (^ we plot the dependence of l^+p/IA"!^ on nit/niH for fixed tan(3 = 2 and 
four different /i scales, {mw, 10, 5, 2.5) CeV in model II and model I, respectively. Decreasing 
the scale fi weakens the dependence of the ratio |yl"'"p/|A~p on mn and the contribution of the 
charged Higgs bosons to gets small. The lower limit of the Higgs mass is sensitive 

to the scale fi and it increases with decreasing in model II, (Table (0)). However, the Higgs 
mass has no lower bound in model I. 

Fig. (^ and (|^) show the dependence of lA+p/IA"!^ on tan/3 for fixed mu = 500 GeV. 
This ratio is sensitive only to small tan/S values. The fi scale regulates the lower limit of tanj3 for 
model I and model II in an opposite way. Further, the effect of the charged Higgs contribution 
becomes weak for large tanf3 values. 

In fig. (|^) and (0) we present the /i scale dependence of |y4+p/|y4^p for the SM and 2HDM 
with tanjS = 10 and two different mass values mn = 500 CeV, 800 CeV for model II and model 
I, respectively. We find, that for model II the smaller the value of mn, the less dependent is 
the ratio on /i. Model I does not allow us to discriminate between the SM or different values 
of niH used as expected (see e.g. fig. (5)). 



set 1 


set 2 


set 3 


As = 370 MeV 
rrib = 5.03 GeV 


A, = 480 MeV 
TTib = 4.91 GeV 


A, = 590 MeV 
nib = 4.79 GeV 



Table 3: The parameter sets of the bound state model, (mf,, A^). 
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The lowest order result of in is obtained by setting /i = niw and it is 0.30 

in the SM for set 3. It reaches 0.85 at /i = 2.5 GeV. Varying fi in the range 2.5 GeV < fi < 
10.0 GeV, is changing in the range 0.60 < < 0.85, resulting in an 

uncertainty of (|^+|2/|^-|2)|^^5 ^ceV) ~ ±35% in the SM. Now we give an example to compare 
the dependence of the ratio |y4"^p/|y4~p on the scale fi and the 2HDM parameters: In model 
II, for mn = 500 GeV and tan(3 > 2, the lowest order result of the ratio is 0.40 

and it enhances up to 0.50 with decreasing tan(3. However, at /i = 2.5 GeV the ratio reaches 
0.86 and the uncertainty due to the extended Higgs sector is weaker than the one due to the 
scale /i. In model I, the behaviour is the same. 

This shows, that the ratio |yl"^p/|74~p is quite sensitive to QCD corrections and this strong 
/i dependence makes the analysis of the 2HDM parameters mu and tan(3 for the given ex- 
perimental value of the ratio difficult. However, we believe, that the strong /i 
dependence will be reduced with the addition the of next to leading order (NLO) calculation, 
and the analysis on the parameters will be more reliable. Note, that a similar analysis for the 
decay 6 57 is given in 0. 

In addition, there is another uncertainty due to the different parameter sets (Table (^). 
For set 2 |A+|V|A-|2 = 0.36 and for set 1 \A^\'^ /\A'\'^ = 0.40 in the SM and in the lowest 
order of a^. It follows, that the larger rrib (smaller A^), the larger the ratio. We further see that 
having increased mt, the ratio lA+p/IA"^ becomes less sensitive to the scale /i. This ratio 
essentially changes when QCD corrections are taken into account. In the lowest order of a^, 
and A~ depend both on the one particle reducible part (IPR), proportional to O;" = C, and 



in addition A~ contains the one particle irreducible part (IPI), proportional to C2, see eq. ([T0|). 
If we include QCD corrections to the considered ratio, the contribution of C^-^-^ dominates over 
the IPI sector and the value of the ratio increases. This means, that the values of A'^ and A~ 
come close to each other and it can be explained as a cancellation of the IPI sector. 

Now we continue to analyze the Br displayed in a series of figures 8-13. In fig. 8-11, we 
present the ^ and tan(3 dependencies of the Br. Decreasing tuh, the Br increases in model 
II, however, the behaviour is opposite in model I. On the other hand the Br is sensitive to 
small tanjS. For large values of tanjS , the dependence of the Br on tanj3 becomes weak in 
model II. In model I, the 2HDM result almost coincides with the SM one since the charged 
Higgs contribution is proportional to l/(tan/3)^. Similar to the case of |74"'"p/|A~p, the Br is 
strongly dependent on the scale fi, see figs (12-13)). It is enhanced for small values of fi. For 
parameter set 3, the lowest order result is 3.6 ■ 10"'^ in the SM. It increases up to 6.8 ■ lO"'' 
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at /i = 2.5 GeV. Varying in the range 2.5 GeV < jj. < 10.0 GeV, the Br changes between 
5.0 ■ 10~^ < Br < 6.9 ■ 10^^ and this gives an uncertainty ABr/ Br{^ = 5 GeV) ~ ±30% for 
the SM. For set 1 (2), the Br increases up to 1.7 ■ 10^^ (1.0 ■ 10^^) and the uncertainty is also 
increases, as due to the scale /i dependence, namely 39 (35) %. This behaviour of the Br results 
mainly from the l/Ag dependence in amplitudes. 

With the addition of extra Higgs contribution, the uncertainty due to 2HDM parameters 
rriH and tanj3 in the Br becomes large like the one coming from the scale /i. Now we will give 
an example to see the effect of the 2HDM parameters on the Br by choosing set 3. In model 
II for rriH = 500 GeV, the lowest order result of the Br is 5.8 ■ 10^'^ for tan(3 > 2 (see fig. 10) 
and it reaches 1.1 ■ 10~^ for smaller tan(3. For comparison, the value in the SM is 3.6 ■ 10~^. 
kt n = 2.5 GeV the Br reaches 6.9 ■ 10"^ (9.0 ■ 10^^) in the SM (2HDM).. This shows, that 
the Br is also sensitive to the extra Higgs contribution. For mu = 500 GeV and /i = 2.5 GeV , 
the Br in the 2HDM model II is enhanced ~ 30% compared to the SM. In model I, there is a 
suppression due to the extra Higgs contribution compared to the SM (see fig. (11)), however, 
the Br is still sensitive to the scale /i (see fig. (13)) and the 2HDM parameters. 

We complete this section by taking the Oj type long distance effects {LDor) for both the 
ratio lA'^p/IA^p and the Br into account. The LDpr contribution to the CP-odd and CP- 
even y4+ amplitudes has recently been calculated with the help of the Vector Meson Dominance 
model (VMD) [P] and it was shown, that the influence on the amplitudes is destructive. With 
the addition of the LDq^ effects, the amplitudes entering |y4+p/|A~p and the Br are now given 
as 

A~ = ^SD + ^Ido, ' (24) 

where Af^, are the short distance amplitudes we took into account in the previous sections 
(eq. 10). The LDq-j amplitudes A^^j^^ are defined as 

Alno, = V2^^FmUmt^C^7"if^) , (25) 

where /(/,(0) = 0.18 GeV is the decay constant of meson at zero momentum, -Fi(O) is the 
extrapolated Bs (p form factor (for details see 0). 

In fig. ( |T4| - |T9| ) we present the rriH, tan(3 and fi dependencies of the ratio lA+p/IA"!^ 



and the Br with the addition of LDo^ effects for set 3. Here we use -Fi(O) = 0.16 [^]. It 
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can be shown, that the value of the ratio |y4+p/|y4~p decreases with the addition of LDqj 
effects. However, while the scale is decreasing, the effect of the LDo^ contribution on the 
ratio is also decreasing, see fig. 4 and 15. On the other hand, the uncertainty 

resulting from varying fi between 2.5 GeV < fi < 10.0 GeV has increased compared to the case 
without the inclusion of the LDq^j amplitudes: 

±40% . 



A(|A+ 








y\A-\ 




5 GeV) 



The Br decreases with the addition of LDqj effects, since the effect is destructive. The /i scale 
uncertainty of the Br is smaller compared to the case where no LD effect is included and we 
get for the range 2.5 GeV < /z < 10.0 GeV in the SM 

A fir 



Br{ix = 5 GeV) 



±27% . (26) 



The present experimental limit on the decay B^ — > 77 is 

Br{B, 77) < 1.48 ■ 10-^ , (27) 

which is far from the theoretical results. By varying the parameters fi,mH,tanj3, (m^, A^), it 
is possible to enhance the Br up to 2.1 (2.5) ■ 10~^ in model II for mu = 480 GeV and large 
tan(3, where the possible maximal value in the SM is 1.4 (1.7) • 10^®, both at /i = 2.5 GeV and 
for set 1. The numbers in parentheses correspond to the case where no LD07 effects are taken 
into account. 

LLog calculations show that the Br strongly depends on the scale fi. This strong dependence 
will disappear with the addition of NLO QCD corrections. From NLO b s'-f decay, the choice 
of /i = mb/2 in the LLog expression reproduces effectively the NLO result, so one suggest that 
this may work also for the b 577 decay. An additional theoretical uncertainty arises from 
the poor knowledge of the Bs bound state effects. 

We find that the Br increases with the addition of the extra Higgs contribution and even 
at the scale fi = 2.5 GeV this value is ~ 2 orders of magnitude smaller than the present 
experimental upper bound. The other possibility for an enhancement of the Br is the extension 
of the Higgs sector. This forces us to think of further models like MSSM,...etc. and FCNC 
Bs 77 decay will be an efficient tool to search for new physics beyond the SM. 
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Figure 1: tanf3 as a function of the mass ttih for fixed Cy = —0.4049 in the model II of the 
2HDM. Here sohd curves correspond to the scale /i — 2.5 GeV, dashed curves to /j, — 5 GeV 
and small dashed curves to = 10 GeV. 
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Figure 2: mt/niH dependence of the ratio p for set 3 and tanP = 2. Here, solid lines 

(curves) correspond to the SM (model II 2HDM) at fx = mw, long dashed hnes (curves) to SM 
(model II 2HDM) at // = 10 GeV , medium dashed hnes (curves) to SM (model II 2HDM) at 
= 5 GeV and small dashed hnes (curves) to SM (model II 2HDM) at /.( = 2.5 GeV . 
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Figure 3: Same as fig 2, but in model I. 
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Figure 4: tanf3 dependence of the ratio p for set 3 for rriH = 500 GeV. Here, solid 

lines (curves) correspond to the SM (model II 2HDM) at /x = mw, long dashed lines (curves) to 
SM (model II 2HDM) at = 10 GeV, medium dashed hues (curves) to SM (model II 2HDM) 
at // = 5 GeV and small dashed lines (curves) to SM (model II 2HDM) at // = 2.5 GeV. 
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Figure 5: Same as fig 4, but in model I 2HDM. 
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Figure 6: The scale dependence of the ratio for set 3 in the SM and for 2 values 

of niH — 500, 800 GeV with tanP — 10. Here, sohd lines (curves) correspond to the SM at 
II — mw (at arbitrary scale) , long dashed lines (curves) to model II with mn — 800 GeV at 
/X = mw (at arbitrary /i scale), and small dashed lines (curves) model II with mn — 500 GeV 
at II — mw (at arbitrary scale) . 
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Figure 7: Same as fig 6, but for model I. All curves coincide within errors. 
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Figure 8: rrit/mH dependence of the branching ratio Br for set 3 with tanP = 2. Here, sohd 
hnes (curves) correspond to the SM (model II 2HDM) at = mw, long dashed lines (curves) 
to the SM (model II 2HDM) at // = 10 GeV, medium dashed hnes (curves) to the SM (model 
II 2HDM) at /X = 5 GeV and small dashed hnes (curves) to the SM (model II 2HDM) at 
H = 2.5 GeV. 
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Figure 9: Same as fig 8, but in model I 2HDM. 
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Figure 10: tanf3 dependence of the Br for set 3 with itlh = 500 GeV. Here, soHd hnes (curves) 
correspond to the SM (model II 2HDM) at = mw, long dashed lines (curves) to the SM 
(model II 2HDM) at = 10 GeV, medium dashed lines (curves) to the SM (model II 2HDM) 
at // = 5 GeV and small dashed lines (curves) to the SM (model II 2HDM) a,t fj, — 2.5 GeV. 
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Figure 11: Same as fig 10, but in model I 2HDM. 
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Figure 12: The scale dependence of the Br for set 3 and for 2 values of mn = 500, 800 GeV 
with tan(3 — 10. Here, solid lines (curves) correspond to the SM at /i — mw (at arbitrary n 
scale) , dashed hues (curves) to model II 2HDM with mn — 800 GeV a,t /j, — mw (at arbitrary 
jjL scale), and small dashed hues (curves) model II 2HDM with mn — 500 GeV ai /i — mw (at 
arbitrary /i scale). 




Figure 13: Same as fig 12, but for model I 2HDM. All curves coincide within errors. 
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Figure 14: mt/mH dependence of the ratio for set 3 at tarif] = 10 with the addition 

of I/i^Oy effects. Here, solid lines (curves) correspond to the SM (model II 2HDM) at = mw, 
long dashed lines (curves) to the SM (model II 2HDM) at = 10 GeV, medium dashed hues 
(curves) to the SM (model II 2HDM) at = 5 GeV and small dashed hues (curves) to the SM 
(model II 2HDM) at = 2.5 GeV. 
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Figure 15: tanP dependence of the ratio l/l"*"!^/!/!"!^ for set 3 for mn = 500 GeV with the 
addition of LDor effects. Here, sohd hues (curves) correspond to the SM (model II 2HDM) 
ai ijl — mw, long dashed hues (curves) to the SM (model II 2HDM) at = 10 GeV, medium 
dashed hues (curves) to the SM (model 11 2HDM) ji — b GeV and small dashed lines (curves) 
to the SM (model II 2HDM) at = 2.5 GeV. 
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Figure 16: The scale dependence of the ratio for set 3 and for 2 values of mn = 

500, 800 GeV at tan(3 = 10 with the addition of LDq^ effects. Here, solid lines (curves) 
correspond to the SM at /i — mw (at arbitrary /i scale), dashed lines (curves) to model II 
2HDM with uih — 800 GeV a,t pL — mw (at arbitrary scale), and small dashed lines (curves) 
to model II 2HDM with uih — 500 GeV a,t /i — mw (at arbitrary fj, scale). 
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Figure 17: mt/mu dependence of the Br for set 3 at tanP = 10 including LDq^ effects. Here, 
solid hues (curves) correspond to the SM (model II 2HDM) at = mw, long dashed hues 
(curves) to the SM (model II 2HDM) at /.t = 10 GeV, medium dashed hues (curves) to the SM 
(model II 2HDM) at = 5 GeV and small dashed lines (curves) to the SM (model II 2HDM) 
at /X = 2.5 GeV. 
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Figure 18: tan/3 dependence of the Br for set 3 for ttih = 500 GeV including LDqy effects. 
Here, solid lines (curves) correspond to the SM (model II 2HDM) at = m^y, long dashed lines 
(curves) to the SM (model II 2HDM) at = 10 GeV, medium dashed hues (curves) to the SM 
(model II 2HDM) at = 5 GeV and small dashed lines (curves) to the SM (model II 2HDM) 
at = 2.5 GeV. 




Figure 19: The scale dependence of the Br for set 3 and for 2 values of mn = 500, 800 GeV 
at tanl3 — 10 including LDqj effects. Here, sohd hnes (curves) correspond to the SM at 
ji — mw (at arbitrary ji scale) , dashed lines (curves) to model II 2HDM with tuh — 800 GeV 
at /X = mvi/(at arbitrary ^ scale) and small dashed lines (curves) to model II 2HDM with 
rriH — 500 GeV at /i — mw (at arbitrary /i scale). 
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